function res = gammap(nu, N)

Arg = [nu-N+1:nu];
res = pi^(N*(N-1)/4)*prod(gamma(Arg/2));